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Abstract 

The accurate determination of atomic masses and the high-precision measurement of the 
bound-electron g factor are prerequisites for the determination of the electron mass, which 
is one of the fundamental constants of nature. In the 2002 CODATA adjustment [P. J. Mohr 
and B. N. Taylor, Rev. Mod. Phys. 77, 1 (2005)], the values of the electron mass and the 
electron-proton mass ratio are mainly based on g factor measurements in combination with 
atomic mass measurements. In this paper, we briefly discuss the prospects for obtaining 
other fundamental information from bound-electron g factor measurements, we present 
some details of a recent investigation of two-loop binding collections to the g factor, and we 
also investigate the radiative corrections in the limit of highly excited Rydberg S states with 
a long lifetime, where the g factor might be explored using a double resonance experiment. 

PACS Nos.: 31.30.Jv, 12.20.Ds, ll.10.St 
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1 Introduction 

The central equation for the determination of the electron mass m c from g factor measurements 
reads 

w c 9 M 

m e = — m ion , (1) 

u L 2q 

where uj c is the cyclotron frequency of the ion, cu L the Larmor spin precession frequency, q 
the ion charge, and m ion its mass. The quantity e = —\e\ is the elementary charge, and g is 
the bound-electron g factor. In most practical applications, the ion is hydrogenlike, and the 
frequency ratio lu c /lul can be determined very accurately in a Penning trap JUll]. 

Equation dU) may now be interpreted in different ways: 

• The ratio m e /m ion is immediately accessible, provided we assume that quantum electro- 
dynamic theory holds for g. Provided the ratio m ion /m p (with the proton mass m p ) is 
also available to sufficient accuracy, the electron to proton mass ratio m c /m p can be de- 
termined by multiplication m c /m ion x m ion /m p . In the recent CODATA adjustment lH, 
the ratio m e /m p has been determined using two measurements involving 12 C. 



*Dedicated to Professor H.-Jiirgen Kluge on the occasion of the 65 birthday 
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• Let us suppose that m ion is known to sufficient accuracy. Assuming that quantum electro- 
dynamic theory holds for g, we may then determine m e from the measurement JHIHEl. 



• The g factor depends on the reduced mass of the electron-ion two-particle system. An 
accurate measurement of g can therefore yield an independent verification of the isotopic 
nuclear mass difference, provided that the masses of the ions have been determined be- 
forehand to sufficient accuracy [J. 

• Direct access to the electron g factor in a weak external magnetic field depends on the 
property of the nucleus having zero spin. According to a relatively recent proposal 10,01, 
the measurement of a g factor for a nucleus with non-zero spin can be used to infer the 
nuclear g factor, provided the purely electronic part of the g factor is known to sufficient 
accuracy from other measurements. 

• There is also a proposal for measuring g factors in lithiumlike systems, and theoretical 
work in this direction has been undertaken II 1011 . Provided the contribution due to electron- 
electron correlation can be tackled to sufficient accuracy, a measurement of the g factor 
in lithiumlike systems could give access to the nuclear-size effect, which in turn can be 
used as an additional input for other determinations of fundamental constants. 

• Finally, provided the mass m- lon of a high-Z ion is known to sufficient accuracy and m e 
is taken from g factor measurements at lower nuclear charge number, the high-Z ex- 
perimental result for g may be compared to a theoretical prediction, yielding a test of 
quantum electrodynamics for a bound particle subject to an external magnetic field and a 
strong Coulomb field (see, e.g., Sec. 2.2 of [ 11]). Alternatively, one may invert the rela- 
tion g = g(a) to solve for the fine-structure constant (important precondition: knowledge 
of nuclear size effect) JslLLJ]. The feasibility of the latter endeavour in various ranges of 
nuclear charge numbers will be discussed in the current article. 

These examples illustrate the rich physics implied by g factor measurements in combination 
with the determination of atomic masses via Penning traps. Indeed, the g factor is a tremendous 
source of information regarding fundamental constants, fundamental interactions and nuclear 
properties. 

This paper is organized as follows. In Sec.|2j we briefly discuss the importance and the status 
of atomic mass measurements for further advances. In Sec.|3l we describe a few details of two 
recent investigations Q, |g] regarding one- and two-loop binding corrections to the g factor, 
and in Sec.|H we discuss the asymptotics of the corrections for high quantum numbers, with a 
partially surprising result, before dwelling on connections of the g factor to nuclear effects and 
the fine- structure constant in Sec.|5J Conclusions are drawn in Sec.|6j An Appendix is devoted 
to the current status of the free-electron anomaly. 



2 Atomic Mass Measurements - Present and Future 



A review of the current status of atomic mass measurements can be in found in Ref. II 1 311 



Experimental details regarding modern atomic mass measurements, with a special emphasis 
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on hydrogenlike ions, can be found in Refs. 11141 Il5ll . Regarding the current status of mass 
measurements, one may point out that some of the masses of S, Kr and Xe ions have recently 
been determined with an accuracy of better than 1 part in 10 10 (Ref. [ 16]). For molecular ions, 
the accuracy has recently been pushed below 1CT 11 lITvIl . 

Recent measurements for the hydrogenlike ions 24 Mg 11+ and 26 Mg 11+ (Ref. Q) and 40 Ca 19+ 
(Ref. HH), as well as for the lithiumlike ion 40 Ca 17+ (Ref. [18]) have reached an accuracy of 
about 5 x 1CT 10 . These experiments pave the way for accurate determinations of fundamen- 
tal constants using g factor measurements in these systems. At the University of Mainz 
(MATS collaboration) and at the University of Stockholm JH] (SMILE-TRAP), there are plans 
to significantly extend and enhance atomic mass measurements (including many more isotopes 
and nuclei) over the next few years, with accuracies below 1 part in 10 11 or even 10 12 . Eventu- 
ally, one may even hope to determine the nuclear size effect of a specific ion by "weighing" the 
Lamb shift of the ground state. In the same context, one may point out that the masses of differ- 
ent charge states of ions are determined vice versa by adding and subtracting binding energies. 
This implies, e.g., that the mass of 12 C 5+ in terms of the mass of neutral carbon, m( 12 C) = 12 u, 
is given by 

m( 12 C 5+ ) = m( 12 C) -5m e + C - 2 £ B , (2) 

where E B = 579.835(1) x 10~ 9 uc 2 is the cumulative binding energy for all 5 electrons lEoTl . 
This relation has proven useful in the determination of the electron mass 0]. 

In order to make a comparison to the accuracy of the free-electron determination of a, it is per- 
haps useful to remember that in the seminal work lEHl . the free-electron and positron anomaly 
has been determined to an accuracy 4 x 10~ 9 . This translates into a level of accuracy of about 
4 x 1CT 12 for the g factor itself. The accuracy of the current value of a is 4 x 10~ 9 yD. 



3 Calculation of the Bound-Electron g Factor 

The bound-electron g factor measures the energy change of a bound electron (hydrogenlike ion, 
spinless nucleus) under a quantal change in the projection of the total angular momentum with 
respect to an axis defined by a (weak) external magnetic field. In this sense, the g factor of a 
bound electron should rather be termed the gj factor (according to the Lande formulation). 

However, for S states, the total angular momentum number is equal to the spin quantum number, 
and therefore it has been common terminology not to distinguish the notation for g and gj. 

For a general hydrogenic state, the Dirac-theory g factor, denoted gn, reads (see |0] and refer- 
ences therein) 

k ( E nK r 



Here, E n j is the Dirac energy, and the quantum numbers n, j and n have their usual meaning. 
Throughout this article, we use natural units with h = c = e = 1. 

For S, P and D states, Eq. © leads to the following expressions (we here expand the bound- 
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state energy in powers of Za), 



, „ ^ 2(Zaf (Za) 4 (2 1 , 

9D(n5l/2 ) = 2 __L_L_i_^___), (4 a) 

9d( „p 1/2) = | i,, (4b) 

, D , 4 8(Za) 2 (Za) 4 /4 2 , 
3D(nP 3 / 2 ) = - - i - — - — , (4c) 

3 15 rr n 6 \15 on' 

4 8(Za) 2 (Za) 4 ( 4 2 . 

<to n£> 3 /2 = - - -±-4r - - - — , (4d) 

5 15 n z rr \ 15 on ,' 



, _ . 6 18 (Za) 2 (Za) 4 / 6 27 » 
gv{nD h/2 ) = - - \ J - — - — . (4e) 

5 35 rr n J V 35 70 n ' 



The above formulas illustrate the in principle well-known fact that the bound-electron g factor 
would be different from the free-electron Dirac value g = 2, even for S states and even in the 
absence of quantum electrodynamic loop corrections. 

We now briefly summarize the results of recent investigations @] of the bound-electron g 
factor, which is based on nonrelativistic quantum electrodynamics (NRQED). The central result 
of this investigation is the following semi-analytic expansion in powers of Za and ln(Za) for 
the bound-electron g factor (nS state) in the non-recoil and pointlike-nucleus limit (for recoil 
effects see e.g. Ref. [ 22]): 



Breit (1928), Dirac theory 



+ 1 1 2 x \{ i + ^&Pi + ^-^ 4i in ^)~ 2 ] +a4 ^ + 



one-loop correction 



+ (^) 2 {-0.656958 (l + ^f^J + M(Za)~ 2 ] + &4o| + 0(Za) 5 | 
^ .. ' 



two-loop correction 



+ 0{a z ) . (5) 

This expansion is valid through the order of two loops (terms of order a 3 are neglected). The 
notation is in part inspired by the usual conventions for Lamb-shift coefficients 12311 : the (lower 
case) a terms denote the one-loop effects, with a^j denoting the coefficient of a term propor- 
tional to a {Za) k \r^[{Za)^ 2 }. The b terms denote the two-loop corrections, with bkj multiply- 
ing a term proportional to a 2 (Za) h ln J, [(Za)~ 2 ]. In JUI^], complete results are derived for the 
coefficients 041, 040, 641 and 640, valid for arbitrary excited S states in hydrogenlike systems. 
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Figure 1: One-loop, two-vertex scattering diagrams that correspond to the one- 
loop part of the effective operators Eq. (B16) and (B17) of Ref. [6]. The zigzag line 
denotes the interaction with the external field, whereas the dashed lines denote the 
Coulomb photons. The two-loop part of these effective operators is generated by 
diagrams with one more virtual photon, with electron-photon vertices to be inserted 
at all topologically distinguishable positions in the electron lines of diagrams (a), 
(b) and (c). 



In Eq. ©, the term underlined by "Breit (1928), Dirac theory" corresponds to the prediction 
of relativistic atomic theory, including the relativistic corrections to the wave function H2411 . By 
contrast, the term - (2 x ~) in the expression underlined by "one-loop correction" gives just the 
leading (Schwinger) correction to the anomalous magnetic moment of a free electron. This latter 
effect is modified here by additional binding corrections to the one-loop correction, which give 
rise e.g. to terms of order a (Za) 2 and higher (in Za). Perhaps, it is also worth clarifying that 
the term —0.656958 is just twice the two-loop contribution to the anomalous magnetic moment 
of a free electron, which is usually quoted as (^) 2 (—0.328479) in the literature. 

Explicit results for the coefficients in ©, restricted to the one-loop self-energy, read H] 

a 41 (nS) = ^ , (6a) 
y 

73 5 8 8 

aAo(nS) = In ko(nS) In k-i(nS) . (6b) 

ov ; 54 24n 9 V ; 3 V ; 

Here, In ko(nS) is the Bethe logarithm for an nS state, and In ^(nS) is a generalization of the 
Bethe logarithm to a perturbative potential of the form 1/r 3 (see also Table [D below). Vacuum 
polarization adds a further n-independent contribution of (—16/15) to a 40 H25I1 . Higher-order 
binding corrections to the one-loop self-energy contribution to the g factor have been consid- 



ered, e.g., in [26], and for the vacuum-polarization contribution, one may consult, e.g., Ref. H27I1 . 
The results for the two-loop coefficients read 

hi{nS) = ^, (7a) 

b 40 (nS) = ^^--\nk --\nUnS) + — tt 2 - — vr 2 ln2 + — C(3) 

3 V ; 810 90 60 SV ; 

tt^-tt 2 ln2 _l C (3)^ . ( 7b ) 
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Our result for 6 40 includes contributions from all two-loop effects (see Fig. 21 of B28J for the 
diagrams) up to the order a 2 (Za) 4 . The logarithmic term 641 is, however, exclusively related to 
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the two-loop self-energy diagrams. An essential contribution to the one- and two-loop effects is 
given by two-Coulomb-vertex scattering amplitudes (see also Fig.[l]). 

Dependence of In k 3 (nS) on n 
5.5-1 1 

5- 

4.5- 

lnfc 3 

4- 
3.5- 



0123456789 
Principal Quantum Number n 

Figure 2: A plot of the generalized Bethe logarithms In k 3 (nS) as a function of 
the principal quantum number n illustrates the monotonic increase with n. For the 
numerical values, see Table[T] 



4 Asymptotics for High Quantum Numbers 

It is interesting to study the limit of the coefficients 040 and 640 in the limit of highly excited 
states, n — > 00. For the Bethe logarithm In k , such a study has recently been completed (see 
Refs. lEil lioll ). The asymptotics of the generalized Bethe logarithm In k s have not yet been 
determined. We here supplement the numerical result for 8S. In Eq. (72) of Ref. @], results 
have been communicated for S states with n < 7. 



Dependence oi*lnA^(?!5) on l j n 



111 /,".; 




Inverse Principal Quantum Number 1 / n 



Figure 3: A plot of the generalized Bethe logarithms \nks(nS) as a function 1/n 
instead of n indicates consistency with an asymptotic limit lim^oo In k^{nS) = 
10 ±2. 
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Table 1: A table of generalized Bethe logarithms In ks(nS) 
for excited S states. This quantity enters into Eqs. (I^El) 
and (TTbh and characterize the one-loop binding correction 
to the g factor of order a [Zaf and the two-loop correction 
of order a 2 (Za) 4 . All decimals shown are significant. 
n In k^(nS) 

1 3.272 806 545 

2 3.546 018 666 

3 3.881 960 979 

4 4.178 190 961 

5 4.433 243 558 

6 4.654 608 237 

7 4.849 173 615 

8 5.022 275 220 



The result for n = 8 confirms the trend of a monotonic increase of In k 3 with n (see Fig. O. On 
the other hand, based on the general experience regarding the structure of radiative corrections 
in the limit n — > oo, we would expect a constant limit of \nks(nS) for n — > oo. Using an 
extrapolation scheme similar to the one employed in 113 ill , we conjecture the following limit 
(see Fig. 

lim lnk 3 (nS) = 10 ± 2 , (8) 

n— »oo 

It would be very interesting to verify this limit by an explicit calculation, e.g., using the tech- 
niques outlined in Ref. W29I1 . 

Highly excited Rydberg states are characterized by a long lifetime. In a Penning trap, however, 
the confining electric fields would tend to quench transitions to lower- lying levels. One might 
attempt a measurement of a g factor of a Rydberg state via a double-resonance approach, with 
one laser driving the spin flip (Larmor precession frequency) and another being tuned to a 
transition between Rydberg states 113211 . 



5 Bound-Electron g Factor, Nuclear Effects and the Fine- 
Structure Constant 

In Figs.|4]and|5l we indicate three primary sources of the theoretical uncertainty of the bound- 
electron g factor across the entire range of nuclear charge numbers (these are the fine- structure 
constant, higher-order unknown two-loop effects and the nuclear radius). For a determination of 
the fine- structure constant using the bound-electron g factor, the experimental accuracy would 
have to be improved to a value below the corresponding uncertainty curve in Figs. |4| and 
Such a determination would constitute a very important and attractive additional pathway, using 
bound-state quantum electrodynamics, as an alternative to the "usual" determination based on 
the free-electron g factor. 
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Before we dwell further on the fine-structure constant, we briefly discuss the nuclear polariz- 
ability correction to the g factor (see Ref. [33] and Appendix© which represents an additional 
obstacle in the determination of the fine- structure constant from g factor measurements. One 
might hope that it can be accurately understood one day in terms of nuclear models. In Ap- 
pendix |X[ we present an additional nuclear effect (a magnetic susceptibility correction) which 
may also have to be taken into account in an accurate description of the nuclear contributions 
to the bound-electron g factor, especially in the range of medium nuclear charge numbers. 

The further shift/uncertainty of the g factor, caused by the nuclear finite-size effect (nuclear 
volume effect), is typically smaller than the uncertainty of the theoretical prediction for the g 
factor due to higher-order quantum electrodynamic two-loop binding corrections (see Figs. HI 
and In evaluating the uncertainty due to the nuclear radius, we have used the most recent 
values for the root-mean- square (rms) nuclear radii [34]. 



Theoretical Uncertainty of the Bound-Electron g Factor 



Sg 




uncertainty due to a 



two-loop uncertainty 

uncertainty due to nuclear radius 



Nuclear Charge Number Z 



Figure 4: Various sources of theoretical uncertainty for the bound-electron g factor, 
over the entire Z range from hydrogen to uranium. 



In order to investigate the sensitivity of the bound-electron g factor to the fine-structure constant, 
we approximate the g factor by the first two terms in the Za-expansion of the Dirac theory and 
the one-loop correction, and obtain 



Sg 



(Za) 



5a 



(9) 



For a determination of a, it is desirable, in principle, to tune the parameters so that the modulus 
\8g\ for given 5a becomes as large as possible. 

For nuclear charge numbers in the (fictitious) range 5 < Z < 6, the sensitivity of g on a suffers 
from a cancellation of the one-loop against the Dirac binding corrections (see also Figs.|4]and|5l), 
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and we have 



5j_ 

5a 



for Z«5.r. 



(10) 



It would be rather difficult to determine a via a measurement of the g factor in the indicated 
range of nuclear charge numbers. 



Theoretical Uncertainty of the g Factor (Small Z) 



&g ltr 11 




uncertainty due to a 



two-loop uncertainty 

uncertainty due to nuclear radius 



4 5 6 7 
Nuclear Charge Number Z 
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Figure 5: A close-up of Fig.H|in the range of small quantum numbers n illustrates 
that an alternative determination of the fine-structure constant, based on the cur- 
rent theoretical status, would be possible for ionized helium ( 4 He + ) and beryllium 
( 10 Be + ). The 6,7 Li nuclei are not spinless. 



For large Z, one may find a crude approximation to Eq. © by the relation 

4 72 \ S 9\ _ 2 2 5a 

-Z a =>- - [Z a) — . (11) 



5a 3 g 3 a 

The enhancement of the theoretical uncertainty of g with Z is manifest in Fig. |4j In principle, 
one might assume that a measurement at high Z could be more favourable for a determination of 
a than a corresponding experiment in a \ow-Z system. However, as shown in Fig.|H the nuclear 
structure alone currently entails an uncertainty of g that is larger than the uncertainty due to the 
fine- structure constant, for large Z. Also, the uncertainty due to higher-order unknown two-loop 
binding corrections currently represents an obstacle for an alternative determination of a from 
a g factor measurement at high Z. 

Reversing the argument, one may point out that, provided the two-loop uncertainty of the the- 
oretical prediction for large Z can be reduced substantially, one may infer the nuclear radius 
from the measurement of the g factor. Again, going one step further and assuming that the nu- 
clear radius is accurately known from other measurements, e.g., Lamb shift experiments or g 
factor measurements in lithiumlike systems, one may eventually hope to infer the fine- structure 
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constant from a high-Z measurement. This endeavour can thus be interpreted as a rather diffi- 
cult combined effort of theory and experiment, with results not to be expected in the immediate 
future, but providing a very interesting perspective in the medium and long term. In particular, 
this endeavour would depend on a successful evaluation, nonperturbative in Za, of all two-loop 
binding corrections to the bound-electron g factor. 

As Fig. 13 shows, the determination of a based on the bound-electron g factor currently appears 
much more promising for extremely light systems, such as 4 He + . The measurement of g factor, 
however, would definitely have to be carried out with an accuracy better than 1CT 11 in order to 
match the current accuracy for a. Alternatively (see Fig.|U), the planned g factor measurement 
in 40 Ca 19+ could potentially lead to a value of a that matches the accuracy of the free-electron 
value, provided the two-loop uncertainty (higher-order binding corrections) can be reduced and 
provided the accuracy of the atomic mass determination can be enhanced beyond 1CT 10 . With 
current theory, the accuracy of the determination of a from the 40 Ca 19+ measurements is limited 
to an accuracy of about two orders of magnitude less than the free-electron value. 

A final word on the electron mass: For a speculative alternative determination of a in a high-Z 
experiment, the current accuracy of m e , based on the carbon and oxygen measurements |ll|,|2l 
flpilii Eill is sufficient. We recall the values (evaluated using the most recent theory |0]) 

m e ( 12 C 5+ ) =0.000 548 579 909 32(29) u , (12) 
m e ( 16 7+ ) =0.000 548 579 909 60(41) u . (13) 

However, if an alternative determination of a via low-Z measurements is pursued in earnest, 
then it becomes necessary to improve the value of m c beyond the 10~ n threshold. 



6 Conclusions 

In Sec. El we emphasize the importance of current high-precision and upcoming ultra-high 
precision atomic mass measurements for the determination of fundamental physical constants, 
in combination with bound-electron g factor measurements in hydrogenlike systems. One of 
the celebrated achievements connected to g factor measurements lies in the improvement of 
the accuracy of the electron mass by a factor of 4, as compared to the previous value based on 



measurements involving protons and electrons in Penning traps [35]. 

The expansion of the bound-electron g factor in terms of the two most important parameters in 
the non-recoil limit is discussed in Sec.|3] These are the loop expansion parameter a (the fine- 
structure constant) and the Coulomb binding parameter Za, where Z is the nuclear charge num- 
ber. Furthermore, in Sec. we analyze generalized Bethe logarithms, termed In k 3 , which are 
relevant for binding corrections to the g factor, in the limit of large principal quantum number 
(i.e., for highly excited Rydberg states). The calculation of the result In k 3 (8S) = 5.022 275 220 
(see Tabled, facilitates the analysis of the asymptotic limit. The discussion is accompanied by a 



tentative proposal B32H for a double-resonance experiment, to probe the bound-electron g factor 
for highly excited Rydberg states with a long lifetime. In Sec. |5J we discuss prospects for de- 
terminations of nuclear properties, and of the fine-structure constant, based on measurements in 
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various ranges of the nuclear charge number. An alternative measurement of the fine- structure 
constant, of comparable accuracy to the free-electron value, could be accomplished via mea- 
surements at low Z, provided the experimental accuracy of the g factor can be pushed beyond 1 
part in 10 11 , and provided the electron mass can be determined to sufficient accuracy (see also 
Figs. IU and 13). A priori, combined ultra-high precision measurements in 4 He + and 10 Be 3+ ap- 
pear to provide for a viable approach, provided the atomic mass measurements of 4 He and 10 Be 
can reach comparable accuracy (now, the experimental accuracy stands at 1.5 parts in 10 11 for 
4 He, see Ref. [ 13]). The two measurements in He and Be could provide input data for a coupled 
system of equations, to be solved for a and m e . 

By contrast, considerable further theoretical and experimental progress (concerning, e.g., nu- 
clear radii) is required before any such endeavour could be realized in the domain of high nu- 
clear charges. The prerequisites are outlined in Sec. 13 We conclude that even in the absence of 
this progress, prospective measurements at higher Z will yield a rather interesting verification 
of quantum electrodynamics in the high-field domain. 
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A Nuclear Magnetic Susceptibility Correction 

It is well recognized that the nuclear polarizability can shift atomic energy levels or electronic g 
factors 113 311 . Less well known is the influence of nuclear magnetic susceptibility (3m, which can 
be significant for large Z-nuclei. The effective interaction Hamiltonian which defines (3 m is 

5H=^fB\ (14) 

where B is the magnetic field at the nucleus. Here, we estimate (3m on the basis of simple 
assumptions. In particular, we assume that nucleus is a bound system of nonrelativistic nucleons. 
Therefore, the Hamiltonian in the external magnetic field is 

« = E - "p »>*) + E - 3 > § ) < 15 > 

where Hj is the interaction Hamiltonian, which we assume to be lEindepedent. The proton and 
neutron masses are denoted by m p and m n , respectively. The term linear in B gives the nuclear 
magnetic moment. The quadratic term from -k 2 = {p + e A) 2 , 

2 m n 8 m„ 

i— i p i— 1 p 
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gives the magnetic suseptibility (3 M (we denote by e the electron charge, e = — |e|). If we 
assume that the quadrupole moment vanishes or is negligible, then 



P 2 9 JL Zp 2 
4m D 3 6m D 

where (r 2 ) is the mean square charge radius of the nucleus. 

We now turn to the correction to the g factor. The magnetic field in Eq. (fT4l is a sum of an 
external magnetic field B cxt and the field B c \ produced by the bound electrons. This leads to an 
effective additional interaction of the electronic magnetic moment with the external magnetic 
field, 

8H = P M -Bext ■ -Bel = "And ' B c \ (18) 

where jl md is the induced nuclear magnetic moment. In the nonrelativistic limit, the magnetic 
dipole interaction of the nuclear magnetic moment fij with the electron in the S-state results in 
a hyperfine structure Hamiltonian 

5H Ms = --^fl I -a e (5 3 (r c )). (19) 
3 m e 

It is now straightforward to write down the analogous Hamiltonian 5H which describes the 
interaction with the induced nuclear magnetic moment /i in d, 

6H = --^- /2 ind • 3 C (6 3 (r c )) = ^— (3 M 5 cxt • a c (6 3 (r c )) = 8g ( -^-) ^ • B cxt (20) 
Sm c Sm e \2m e J 2 

where 

Sg = ~- (5 3 (r e )) Pm = m 3 e (3 M , (21) 

6 6 Tin 6 

and Pm is defined in Eq. (fTTt . As an example, one may consider 40 Ca with Z = 20 and a 



radius of a/ (r 2 ) = 3.4764(10) fm Il34ll . Using values for the fundamental constants as given 
in 02D, one obtains an estimate for the nuclear susceptibility of (3m = 1-35 • 10 _8 m~ 3 and 
^fitea ~ —1.78 • 10~ n for n = 1, which is much less than the uncertainty due to higher order 
two-loop corrections but important for an accurate understanding of nuclear contributions to the 
bound-electron g factor. 
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